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Abstract 

The rapidity dependence of elliptic, V 2 , and triangular, V 3 , flow coefficients in proton-nucleus (p-i-A) collisions is pre¬ 
dicted in hydrodynamics and in a multi-phase transport model (AMPT). We find that v„ {n = 2,3) on a nucleus side 
is significantly larger than on a proton side and the ratio between the two, weekly depends on the transverse 

momentum of produced particles. 


1. Introduction 

Significant second and third harmonics have been ob¬ 
served in the long-range azimuthal correlations of par¬ 
ticles emitted in ultra-relativistic p-tPb collisions at the 
LHC [1-3] and d-i-Au collisions at RHIC [4]. The results 
can be interpreted as due to collective flow of particles 
in the framework of hydrodynamic models [5-11] or in 
the cascade AMPT model* [13-15]. A different approach 
connects the observed particle correlations with saturation 
effects in the initial state of the collision [16-19]. Observ¬ 
ables related to rapidity dependence of the bulk quanti¬ 
ties can be used to disentangle between the two mech¬ 
anisms. An example of such observable is the average 
transverse momentum of produced particles as a function 
of (pseudo)rapidity. The average transverse momentum is 
predicted to be larger on the Pb-going than on the p-going 
side in the hydrodynamic model [ 20 ], while the reverse 
is expected in the color glass condensate [21] (CGC) ap¬ 
proach. The transverse size of the fireball is larger (and 
it lives longer) on the Pb-side which results in a stronger 
collective flow during the evolution, leading to not only a 
larger transverse flow (and larger {p±}), but also a stronger 
elliptic and triangular flow. Such an asymmetry of inte¬ 
grated elliptic flow in p- and Pb-going sides has been ob¬ 
served by the CMS collaboration [22]. In this letter we 
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*As shown in Ref. [12], the AMPT model generates the signal 
mostly due to the escape mechanism which presumably differs from 
hydrodynamics. 


present a calculation of the relative strength of the elliptic 
and the triangular flow in the Pb-going and p-going sides 
as a function of the transverse momentum in the 3-1-1 di¬ 
mensional ( 3 - 1 -ID) viscous hydrodynamic model and in 
the AMPT model. In both cases we observe a significant 
increase of the elliptic and the triangular flow coefficients 
for rapidities corresponding to the Pb-nucleus direction. 
Our predictions extend to larger rapidities than measured 
by the CMS collaboration, and correspond to acceptance 
of the ALICE muon spectrometer. 

In the next section we present our main results. In sec¬ 
tion 3 we offer some comments and we conclude the paper 
in section 4. 

2. Results 

In this Section we present the elliptic and triangular 
flow coefficients in the proton- and the nucleus-going di¬ 
rections calculated in the 3-i-lD hydrodynamics and the 
AMPT model. 

2.1. Hydrodynamics 

The initial density for the hydrodynamic evolution is 
calculated in the Glauber Monte Carlo model. The en¬ 
tropy is deposited at the nucleon-nucleon collision points 
with a Gaussian profile in the transverse plane [ 6 ]. The 
3 -I- ID hydrodynamic calculations are performed event- 
by-event, with shear viscosity 77/5 = 0.08 and bulk vis¬ 
cosity rj! s = 0.04 for T < 170 MeV. At the freeze-out 
temperature of 150 MeV, particles are emitted statistically 
[23], including corrections due to bulk and shear viscosity 


Preprint submitted to Elsevier 


March 13, 2015 







Figure 1: The elliptic flow coefficient on the proton- (-4 < rj < -2.5) 
(dashed line) and the nucleus-going (2.5 < rj < 4) (solid line) sides 
in 0 - 20% p-)-Pb collisions at yfs - 5.02 TeV, as a function of the 
transverse momentum, from 3-1-ID hydrodynamics. 



Figure 3: The ratio of the Pb-going to p-going values of the ellip¬ 
tic (solid line) and triangular (dashed line) flow coefficients, from the 
3-l-lD hydrodynamic calculation. 

(Figs. 1 and 2). Since the average transverse momen¬ 
tum is expected to be larger on the nucleus-going side 
[ 20 ], the difference for integrated flow coefficients is pre¬ 
dicted to be even larger. The origin of the effect in the 
hydrodynamic model can be linked to the longer lifetime 
of the fireball on the nucleus going side, which results in 
a stronger built up of the collective flow. Moreover, for 
rapidities where the freeze-out happens earlier we expect 
much stronger viscosity correction at freeze-out, reducing 
the flow coefficients [27]. As shown in Fig. 3, the ratio of 
the flow coefficients calculated in the forward and back¬ 
ward rapidity bins weakly depends on p±, in the range 
where the hydrodynamic model applies. 


Figure 2: Same as figure 1 but for the triangular flow. 

in the Cooper-Frye formula [24]. The hydrodynamic sim¬ 
ulations reproduce fairly well the measured elliptic and 
triangular flow [ 6 ] , the mass hierarchy of the elliptic flow 
coefficient and of the average transverse momentum of 
identified parficles [25], and fhe inferferomefry radii [26]. 

The cenfralify 0-20% is defined as evenfs wifh fhe num¬ 
ber of wounded nucleons N^, >13. Charged parficles are 
analysed in fhree bins, fhe forward (Pb-going side) 2.5 < 
77 < 4 and backward (p-going side) -4 < 77 < -2.5, and 
fhe cenfral bin I 77 I < 1. The cenfral bin defines fhe refer¬ 
ence evenf-plane for charged parficles wifh 0.25 < 77 ± < 5 
GeV. The flow coefficienfs Vn{2}{p±) for charged parficles 
in fhe forward and backward bins are calculafed wifh re- 
specf fo fhe reference parficles from fhe cenfral bin. 

The ellipfic and friangular flow as a funcfion of fhe 
fransverse momenfum is larger on fhe nucleus-going side 


2.2. A multi-phase transport model 

The AMPT model wifh fhe siring melling mechanism 
proved lo be very effective in describing various feafures 
of p+Pb, d-i-Au and fhe high-mulliplicily p-i-p inleraclions 
dala [13-15].^ The model is initialized wifh sofl strings 
(soft particles) and minijets (hard particles) from HIJING 
[28]. In the string melting scenario both strings and mini¬ 
jets are converted into quarks and anti-quarks that sub¬ 
sequently undergo elastic scatterings with a given cross- 
section, cr, which is a free parameter.^ It was found that a 
cross-section of 1.5 - 3 mb is sufficient to reproduce the 
data in p-i-p and p-i-Pb collisions at the LHC [13, 14], and 
d-i-Au interactions at RHIC [15]. In this paper we choose 
0-^3 mb. 


^We note that approximately 1-2 elastic collisions per partons 
suffice to describe the p-l-Pb data [13, 14]. 

^The AMPT model with cr = 0 is equivalent to HIJING (plus 
hadronic transport [29]). 
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Figure 4: The AMPT model (with string melting) results for the elliptic flow coefficients on a proton (-4 < ?/ < -2.5) and a nucleus (2.5 < 77 < 4) 
sides, Vj and v^'’, for various centrality classes in p+Pb interactions at = 5.02 TeV, as a function of the transverse momentum, In this plot 
jets contribute to V 2 for higher values of p^. 
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In Fig. 4 we present V 2 in the proton- (-4 < 77 < -2.5) 
and the nucleus-going (2.5 < 77 < 4) directions as a func¬ 
tion of the transverse momentum, px- We performed our 
calculations in two different ways. In the first method 
we calculated three two-particle correlation functions be¬ 
tween two bins (i) -4 < 77 < -2.5 and I 77 I < 1, (ii) 
2.5 < 77 < 4 and I 77 I < 1, and (hi) -4 < 77 < -2.5 and 
2.5 < 77 < 4.^ In this way we can extract 
and V 2 V 2 '’ what allows to calculate and separately. 

In the second method we extract and V 2 ’’v™‘^ as 

above however, in this case we extract calculating V 2 
in I 77 I < 1.2 with the rapidity gap between particles be¬ 
ing two units of rapidity. This allows to extract with 
a good approximation and in fact we checked that both 
methods lead to practically indistinguishable results. In 
the following we show the results obtained using the lat¬ 
ter method. 

As seen in Fig. 4, v^’’ is larger that (except periph¬ 
eral collisions) for all calculated values of px- In periph¬ 
eral collisions the single particle 77 distribution, dNjdr] be¬ 
comes almost symmetric in 77 leading to » v^. 

It is important to see how jets influence both and v^’’. 
We subtracted jets by randomizing the azimuthal angles 
between produced jets. In this case jets do not contribute 


'*In all cases we have large enough rapidity separation between bins 
allowing to neglect correlations from jet cones, etc. 


to the two-particle correlation function at A 0 = tt and 
consequently do not contribute to the extracted values of 
V 2 (provided we have large enough rapidity gap between 
bins). As seen in Fig. 5 this procedure modifies V 2 at 
large transverse momenta however, all qualitative features 
remain unchanged.^ 

In Fig. 6 we show the ratio between v^’’ and as a 
function of transverse momentum calculated in the AMPT 
and AMPT JS (jets subtracted) models. 

Finally we performed our calculations for V3. In Fig. 7 
we show Vj’’ and as a function of transverse momentum 
in 0 - 20% p-i-Pb collisions dl ^[s - 5.02 TeV calculated 
in the AMPT model. In Fig. 8 we present the ratios of 
and for 0-20% centrality class in the AMPT 
model, where a larger ratio for V 3 than V 2 is seen.^ 

3. Comments 

Several comments are in order. 


^We also checked our results using a different method of jet sub¬ 
traction. We calculated the two-particles correlation functions for 
cr = 3 and cr = 0 mb. The latter has no contribution from collective 
physics. Finally we calculated ^v^Cimby — V 2 ( 0 mb )2 that is sensi¬ 
tive to collective physics only. We found that both methods lead to 
practically the same results. 

^We checked that as expected the AMPT JS model (jets subtracted) 
gives almost identical results. 
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Figure 5: Same as Fig. 4 except jets are subtracted (JS) from the two-particle azimuthal correlation function reducing both and Vj*’ for higher 
values of p^. 
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Figure 6: The AMPT and AMPT JS (jets subtracted) results for the ratio as a function of the transverse momentum p^ 
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4. Conclusions 
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Figure 7: Same as figure 2 but from the AMPT model. 
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Figure 8: Same as figure 3 but from the AMPT model. 


We repeated our calculations using different definitions 
of centrality classes. We checked several possibilities in¬ 
cluding: cuts of the multiplicity distributions in I 77 I < 1 , 
2.8 < 77 < 5.1 and 5 < rj < 6, and cuts in the number 
of wounded nucleons, A^part- As expected, we found some 
quantitative differences (on fhe level of 20 % for 0 - 20 % 
cenfralify class) however, all qualifafive fealures remained 
unchanged. 

If would be inferesfing fo perform analogous calcula- 
fions in fhe color glass condensafe framework. V 2 on a 
nucleus side is driven by large x parfons in a nucleus and 
small X parfons in a proton wifh fhe opposife sifuafion for 
V 2 on a profon side. Consequenfly in CGC we expecf a 
nonfrivial dependence of V 2 on rapidify in p-i-Pb collisions 
and if is plausible fhaf V 2 (t]) could serve as fhe decisive 
fesf of fhe initial vs. fhe final slate effecls. 


In conclusion, we predicfed fhe rapidity dependence of 
elliptic and Iriangular flow coefficienls in p-i-Pb collisions 
al fhe LHC energy using fhe AMPT and 3-i-lD hydrody¬ 
namics models. We found lhal bolh V 2 and V 3 in cen- 
fral collisions are significanlly larger on a nucleus side 
(2.5 < 77 < 4) fhan on a proton side (-4 < 77 < -2.5) 
and fhe ratio belween fhe Iwo, v^^(p±)/vn(p±), weakly de¬ 
pends on fhe Iransverse momenlum of produced particles. 
The signal is somehow larger in hydrodynamics fhan in 
fhe AMPT model. We also predicfed fhe cenfralify depen¬ 
dence of fhe effecl and found fhaf already for 40 - 60% 
cenfralify class fhe rafio is consisfenl wifh unify. If was 
furlher observed fhaf fhe rafio weakly depends on various 
mefhods of cenfralify definition in p-i-Pb (for 0 - 20% cen- 
Iralily class). Finally, we performed our calculalions wifh 
and wifhouf jel confribufion (by randomizing azimulhal 
angle belween produced jels in AMPT) and found very 
lilfle effecl on fhe rafio whereas fhe individual V 2 coef¬ 
ficienls are obviously slrongly modified al larger pj^. If 
would be inferesfing fo perform analogous calculalions in 
fhe initial slale models of p-i-A inleraclions, where a non¬ 
frivial V 2 dependence on (pseudo)rapidity is expecled. We 
hope our resulls will provide a slronger fesf of fhe collec¬ 
tive dynamics in p-i-A collisions. 
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